External Shock Model for Gamma-Ray 
Bursts during the Prompt Phase 



Charles D. Dermer^ 

Naval Research Laboratory, Code 7653, Washington, DC 20375-5352 



Abstract. The hard X-ray and 7-ray phenomenology of gamma-ray bursts (GRBs) 
can be explained by an external shock model where a single relativistic blast wave 
interacts with the surrounding medium. Besides reproducing the generic spectral be- 
havior of GRB profiles, the external shock model provides quantitative fits to the peak 
flux distribution, the > f s t^Q duration distribution, and the distribution of the peaks 
Epk of the iyF,y spectra of GRBs measured with BATSE. The apparent paradox be- 
tween a relativistic beaming scenario and the empirical finding that Epk values are 
preferentially measured within the triggering range of a GRB detector is resolved by 
this model when blast wave physics and detector triggering criteria are taken into ac- 
count. Some surprising implications follow, namely that the fireball event rate is ~ 1 
per 10** years per Milky Way galaxy for unbeamed sources, and proportionally more if 
fireball outflows are coUimated. This is ~ 3 orders of magnitude larger than normally 
quoted. Most of the clean and dirty fireball transients are undetected due to telescope 
sensitivity and design limitations. 

Strongly variable GRB time histories with good radiative efficiencies are possible 
because of the strongly enhanced emissions when a blast wave interacts with den- 
sity inhomogeneities located nearly along the line-of-sight to the observer. Arguments 
against short timescale variability in an external shock model are answered, and difh- 
culties in an internal shock/colliding shell model are mentioned. 



I INTRODUCTION 

An important question in GRB studies is whether the GRB engine produces a 
single impulsive collapse and ejection event, or instead operates over a period of 
time much longer than the ~ ms dynamical time scale of the central engine. In the 
external shock model [1,2], a single relativistic shell is ejected by the GRB engine 
and energized by interactions with the surrounding medium. Variability in the light 
curves is attributed to interactions with an inhomogeneous surrounding medium. 
In the colliding shell (or internal shock) model [3,4], collisions between a succession 
of shells in a relativistic wind are thought to produce the variability observed in 
GRB light curves. If a conclusive resolution to this problem is obtained, then 
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physical information can be extracted directly from GRB light curves. In the case 
of the external shock model, variations in GRB light curves reveal the distribution 

of circumstcllar material near the sources of GRBs. In the case of the internal shock 
model, GRB light curves reflect the structure of and accretion processes operating 
within the putative disk of material that is accreted by the newly formed collapsed 
object to energize the relativistic wind. 

Here we review work focusing on the external shock model in the prompt 7- 
ray luminous phase. We find that the extensive phenomenology of GRBs can be 
explained with this model, so that the addition of multiple relativistic shells and 
the numerous parameters associated with a hybrid internal/external shock model 
are unnecessary. The fewer number of free parameters in the external shock model 
places definite constraints on the number and type of fireballs needed to explain 
GRB statistics. The most important implication is that classes of clean and dirty 
fireballs with well-defined properties must exist, and that the fireball event rate is 
much larger than previously estimated on the basis of detected GRBs. 

II NUMERICAL SIMULATION OF LIGHT CURVES 

When a relativistic blast wave with Lorentz factor F encounters an external 
medium, charged particles will be captured by the blast-wave shell even if the shell 
has only a very weak entrained magnetic field. A captured particle in the shell 
frame gets Lorentz factor F. This internal energy derives from the directed energy 
of the relativistic shell, causing the shell to decelerate. 

We have developed a numerical simulation model [5,6] for a GRB blast wave 
that interacts with an external medium. The model treats synchrotron, Compton, 
and adiabatic processes, and blast-wave deceleration is self-consistently calculated. 
The parameters that enter the numerical model are those of the standard blast 
wave model. The macroscopic variables are the implied isotropic energy release 
E = 10^^-E54 ergs/(47r sr) and the initial Lorentz factor Fq = 3OOF300 of the blast 
wave. The environmental variables are the external density n{x) oc n^x^"^, where 
X is the distance from the center of the explosion. We let no = 10^n2 cm~^ and 
consider a uniform surrounding medium (r^ = 0) . (Inhomogeneities in the external 
medium are considered in §V.) We also let the opening half- angle of the outflow 
ip = 10°, corresponding to a beaming factor / = 0.76%. As long as -i/^ ^ Fg^, 
the coUimation has little effect on 7-ray emission during the prompt phase if the 
observer's line-of-sight falls within an angle 9 ^T^^ of the jet axis. 

The microscopic variables are the fraction of energy ee that is transferred from 
the swept-up protons to the swept-up electrons, and the injection index p of the as- 
sumed power-law electron energy distribution. A parameter emax is defined in terms 
of a maximum Lorentz factor 7inax obtained by balancing the minimum acceleration 
time scale and the synchrotron loss time scale, giving 7inax = 4 x 10'''einax/[-B(G)]-'^/^. 
The magnetic field B is specified by a magnetic- field parameter through the re- 
lation 57(87r) = AeBmpC^n{x)(5{T'^ - F), where (5c = {I - T-'^f/^c is the speed 




FIGURE 1. Calculated GRB light curves (thick curves) and photon spectral indices (thin 

curves) due to an external shock interacting with a uniform surrounding medium (left) arc shown 
in the left figure. The 20-50 and 100-300 kcV comit rates of a typical GRB with a smooth profile 
(GRB 990220; BATSE trigger 7403) are shown in the middle and right figures, respectively. 

of the blast wave. Standard values used here are = 0.5, p = 2.5, Cmax = 1, 
and 63 = 10~^. The low value of is required to avoid forming cooling spectra, 
which are not commonly observed in GRBs [7] . We also note that the microscopic 
variables are assumed to be constant in time. 

The left panel in Fig. 1 shows calculations of hght curves and spectral indices 
at different observing energies for a model GRB with standard parameters. For 
comparison, we also show a typical GRB with a smooth light curve. Several effects 
are apparent here. The first is that the generic Fast Rise, Exponential Decay 
(FRED) profile found in some 20-30% of all GRB light curves is reproduced (FRED 
is actually a misnomer, as the decay law is more closely approximated by a power 
law). The second is that the peaks are sharper at higher energies and broader at 
lower energies. Another is a hardness- intensity relation and a hard-to-soft evolution 
of the GRB light curves, so that the well-known correlations are reproduced. A 
prediction of the model is that the peaks are aligned at 7-ray energies, but lag 
at X-ray energies [8]. This prediction seems to be confirmed by observations with 
Beppo-SAX [9] which has spectral coverage in the 2-700 keV range. 

Fig. 2a shows model GRB spectra at different observing times, and Fig. 2b shows 
the calculated relationship between Ep\^, flux, and fiuence. At X-ray energies, the 
photon spectral index approaches a value a ~ 2/3, corresponding to the nonthermal 
synchrotron emissivity spectrum from an uncooled electron distribution with a 
low-cncrgy cutoff. The spectrum turns over and approaches the value 1 -t- (p/2) 
associated with a cooling electron distribution at the highest energies. Fig. 2b 
shows that the qualitative behavior of the E'pk-fluence relationship observed in 
GRBs [10] is reproduced. The spectral aging inferred from the decay of E'pk values 
in smooth GRB hght curves is a natural consequence of the external shock model. 

The external shock model therefore accounts for the best established phenomeno- 
logical correlations of FRED-type GRBs [11,8]. 
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FIGURE 2. Generic behavior of a model GRB from an external shock energized by a uniform 
surrounding medium. The left panel shows the broadband X-ray and 7-ray spectra at different 
observing times. The right panel shows the dependence of -Bpk and flux as a function of fluence. 



Ill STATISTICAL PROPERTIES OF GRBS 



Even if beaming is neglected, seven parameters enter into a blast-wave model 
calculation with an assumed uniform surrounding medium. We carried out a pa- 
rameter study [12] showing that GRB observables are most sensitive to the value of 
the initial Lorentz factor (or baryon-loading parameter) Fq of the explosion. The 
typical duration of a GRB in the prompt phase varies as {E /TqUoY^^ at observing 
energies £ ^ Sq. The quantity Sq = i?pk(t = 0) is the photon energy of the peak of 
the I'Fi, spectrum at early times, and £0 oc qnl/^rQ, where g is a parameter related 
to the magnetic field and Lorentz factor of the lowest energy electrons. The power 
rio at photon energy Sq varies as (FQ^'Q/no)^/^. These relations show that the mean 
duration, peak photon energy, and peak power output of a GRB are most sensitive 
to the value of Fq. 

A central criticism of a relativistic beaming scenario has been to explain the 
apparent paradox between a model involving relativistically beamed outflows, and 
observations showing that Ep]^ is narrowly confined to an energy range near a 
few hundred keV. Brainerd's Compton attenuation model [13], for example, was 
specifically designed to account for this fact, but the large column densities required 
by this model make it unable to explain rapid variability in GRB light curves [14]. 
The beaming paradox is resolved by the external shock model [8] when the spectral 
behavior implied by blast wave physics is convolved with detector response. A 
dirty fireball with Fq <^ 300 will have a uFiy peak at low energies, and will rarely 
be detected because the blast wave energy is radiated over a long period of time 
(oc Fq^^"^); thus its peak power is very weak (Flo cx: Fq^^). The flux in the BATSE 
range is even lower than imphed by this relation because BATSE would be sensitive 
to only the soft, high-energy portion of the spectrum. A clean fireball, by contrast, 
would produce a brief, very luminous GRB, but BATSE would sample the very 
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FIGURE 3. (a) Prompt 100 keV light curves for blast-wave Lorentz factors To = 1000 (solid 
curve), 300 (dashed), and 100 (dotted). Other parameters are given in the text, (b) Dependence 
of iJpk measured at the deceleration time scale t^ec on the 20 keV - 2 MeV fluence, which is 
integrated from the start of the burst to Stdec- 



hard portion of the spectrum below the uF^ peak where the received flux is not so 
great. Fig. 3a iUustrates this behavior. Dirty fireballs would rarely trigger BATSE 
because the flux is so weak in the BATSE triggering range, and clean flreballs with 
Fq ^ 300 would be so brief that the total fluence measured within the BATSE 
window would not be sufficient to trigger it. 

Fig. 3b shows the relationship between E'pk and fluence for a model calculation 
when only the parameter Fq is varied. When Fq 500, the external shock model 
predicts a positive correlation between i?pk and fluence, as has been recently re- 
ported [16]. The dirty flreballs with Fq 100 would not normally be detected and, 
as just described, there would also be biases against detecting the clean flreballs 
with Fq 3> 300. It is necessary, however, to convolve temporal and spectral model 
results through a simulation of the detector response before drawing conclusions 
about the viability of the model. 

The BATSE instrument has provided the largest and most uniform data base on 
GRBs. It nominally triggers on 64, 256, and 1024 ms timescales when the flux in at 
least two detectors exceeds 5.5(7 over background. The data points in Fig. 4 show 
the peak photon-flux size distribution, the t^Q duration and the i?pk distributions 
measured with BATSE. The observable tso is the time interval over which the 
integrated counts range from 25% to 75% of the total counts over background. For 
comparison with statistical data, we developed an analytic model for the temporally 
evolving GRB spectrum [12] based on the detailed numerical calculations. To make 
a valid comparison between the external shock model and the observed statistical 
properties of GRBs, we have modeled detector triggering criteria. Model results 
were integrated over time to determine if the peak 50-300 keV flux exceeded the 
BATSE threshold so that the simulated BATSE detector would be triggered [15]. 
Trigger efficiencies were explicitly taken into account, which is important for GRBs 
with fluxes near threshold. The underlying assumption of our statistical model is 



that the event rate of fireballs follows the star formation history of the universe 
[17]. 

We [15] found that it was not possible to fit simultaneously the size, t^Q and 
£^pk distributions with a monoparametric model. Broad distributions of explosion 
energy E and initial Lorentz factor Tq are needed to fit these distributions. The 
model fits shown in Fig. 4 are based upon power-law distributions of E and Fq, 
where dN/dE oc E'^-^^ for 10^^ < E(ergs) < 10^^ and dN/dVo oc Fq"'^^ for 
Fq < 260. The upper limit to Fq corresponds to a density of uq = 10^ cm~^; the 
analytic model is degenerate in the quantity ^oFq. 

As can be seen from Fig. 4, the model provides reasonable fits to the peak- 
fiux, i?pk and tso distribution of the long-duration ( ^ 1 s) GRBs. The short hard 
GRBs must arise from a separate component. The implied redshift distribution of 
GRBs detected with BATSE is also shown in Fig. 4. We predict that most GRBs 
detected with BATSE lie in the redshift range 0.2 ^ z ^ 1.2, with a tail of GRBs 
extending to high redshifts. The predicted number and distribution of high-z GRBs 
detected with BATSE is quite uncertain, because the star-formation rate at high 
redshifts is poorly known, and the fit depends on the unproven assumption that 
the comoving space density of fireball transients follows the star formation rate. 
Moreover, the distribution of explosion energies is assumed to be described by a 
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FIGURE 4. Data points give the 3B catalog peak 50-300 keV photon- flux (upper left), t^o 
duration (upper right) and i?pk (lower left) distributions of GRBs measured with BATSE [18,19]. 
Dotted histograms give fits from the external shock model [15] with a range of values of total 
explosion energy and baryon-loading. Lower-right panel shows the redshift distribution predicted 
by this set of parameters. 



power-law function with a discrete cutoff. This distribution might instead have a 
tail extending to very high values. 

IV DIRTY AND CLEAN FIREBALLS 

An overall normalization factor for the fireball event rate per unit comoving 
volume is implied by the joint fits to the statistical properties of GRBs shown 
in Fig. 4. If no beaming is assumed, this normalization corresponds to a local 
event rate of = 440 yr~^ Gpc~^, which is equivalent to a local GRB rate of = 90 
Galactic events per Myr. This is a factor ~ 4000 greater than the result of Wijers 
et al. [20], who fit the combined BATSE/PVO peak-fiux distributions only. This 
difference is due to an approach to GRB statistics where we abandon a standard 
candle assumption for the luminosity and rely on blast wave physics and detector 
response properties to determine whether a fireball transient would be detected 
with BATSE. Most crucially, we do not assume that there is a preference in nature 
to make fireballs with a specific energy release E and baryon-loading parameter 
Fq (which would also entail a typical density of the surrounding medium) that 
would produce radiation that would trigger BATSE; any such assumption is highly 
artificial. 

The consequence of this approach is that fireballs with a wide range of energies 
and baryon-loading parameters are formed in nature, the bulk of which are not 
detected and for which we have no evidence except for the limits implied by surveys 
[21,22]. Only a very few nearby fireball transients with low values of E would be 
detected, and fireballs with Fq ^ 10^ would be invisible to BATSE because most 
of the dirty fireball radiation is emitted at X-ray energies and below. The dirty 
fireball transients have longer durations and lower £'pk values than standard GRBs, 
and are difficult to detect because they arc lost in the glow of the luminous diffuse 
X-ray background for wide field-of-view instruments. The X-ray transient events 
discovered with the Beppo-SAX WFC and reported at this meeting [24] might be 
fireball transients with a baryon load that is large enough that such events would 
not normally trigger a burst detector at hard X-ray energies. The number of clean 
fireball transients is not well constrained, but our results show that there must be 
a break or cutoff in the Fo-distribution at high values of Fq. Clean fireballs have 
shorter durations and E^]^ values extending to MeV and GeV energies, and require 
sensitive, wide field-of-view gamma-ray telescopes to be detected [12,23]. 

If there are many more fireball events than implied by direct observations of 
GRBs, then a number of important implications follow: 

• The hypothesis that ultra-high energy cosmic rays are produced by GRBs 
remains viable. This hypothesis has been questioned [25] in light of redshift 
measurements of GRB counterparts that suggest a much lower event rate 
within the GZK radius than formerly thought. 

• The identification of X-ray hot spots in MlOl with GRB remnants [26] ap- 



pears more probable. These associations seemed unlikely given the event rate 
inferred directly from GRB observations. 

• GRB explosions could leave many more observable Galactic remnants such as 
stellar arcs and HI holes, and produce greater biological effects than has been 
estimated [27]. 

V INHOMOGENEOUS EXTERNAL MEDIUM 

Several arguments have been advanced to the effect that an external shock model 
cannot reproduce the short timescale variability observed in GRB light curves. We 

address these point by point. 

1. An external shock model will display short timescale variability only if the 
radiative efficiency is poor. The analytic argument [28] assumes that density in- 
homogeneities (or "clouds") located at an angle 9 ~ F"^ to the hne-of-sight make 
the dominant contribution to variability. Clouds at ^ <^ F"^ actually make much 
stronger contributions to large-amplitude flux variability because of the combined 
effects of Doppler beaming and the much shorter observer timescale over which 
on-axis clouds radiate their emission [2] . 

2. A condition of local spherical symmetry in radiating blast wave produces pulses 
in light curves which spread with time, contrary to the observations [29]. An external 
shock model breaks the condition of local spherical symmetry if clouds with radius 
r <^ -R/Fo are present, as must be assumed to make the short timescale variability. 
Here R is the distance of the cloud from the explosion center. 

3. A decelerating blast wave produces spreading pulses, contrary to the observa- 
tions. Only the portion of the blast wave that interacts with a cloud experiences 
strong deceleration, and its energy is dissipated by the interaction. The rest of 
the blast wave does not undergo significant deceleration until it intercepts another 
cloud, so no spreading from deceleration results. Thus it is not surprising that 
there is no spreading of peaks in GRB 990123 [30], because different portions of 
the blast wave are producing the distinct pulses and peaks in the light curve. 

4. Gaps and precursors are not possible due to the interference between a large 
number of causally disconnected regions. If there are shells of material from winds 
of GRB progenitor stars, as seems likely if GRB sources are associated with the 
collapse of massive stars, then gaps in the light curves can be formed. 

5. A low-density confining medium will produce a low level of emission unless the 
density contrast between the clouds and the confining medium is very large. First, 
it is not necessary to have a confining medium if the massive star progenitor ejects 
material. Even if there is a low-density confining medium, the standard blast wave 
model implies that this residual emission will be radiated in a different energy band 
than the radiation emitted from the blast- wave/cloud interaction. 

Finally, we note difficulties in a coUiding shell scenario. The efficiency for dissi- 
pating internal energy in a relativistic shell is maximized for collisions between a 
shell and a stationary external medium, and is much poorer in collisions between 



relativistic shells. It is simple to get ^ 10% radiative efficiency in the BATSE band 
for an external shock model, but efficiencies ~ 1% are more likely in an internal 

shock model [31], which calls into question the validity of the internal shock model 
for GRB 970508 [32]. A colliding shell scenario must contend with spreading pro- 
files unless pairs of shells collide only once near the burst source, which would mean 
an additional loss of efficiency. GRBs with widely separated pulses generally have 
uFi, peak photon energies within a factor of 2-3 of each other. This is natural for 
an external shock model, where a blast wave with a single Lorentz factor collides 
with different clouds within the Doppler cone, but requires fine-tuning of the speeds 
between pairs of shells in a colliding shell model. 

VI SUMMARY 

The original motivation for an external shock model was that it provided a simple 
explanation for the mean duration of GRBs [33] . This duration roughly corresponds 
to the time scale tdec where the relativistic shell has swept up a sufficient amount 
of matter to cause the shell to decelerate. For a GRB source at redshift z, 

taeo ^ m + z){^ f' s , (1) 

300 

which is comparable to the mean duration of GRBs observed with BATSE (see 
Fig. 4). This equation does not explain, however, why Fq ~ 300. Wc now know 
the answer to this problem - fireballs do not have to have Fq ~ 300. But if the 
baryon-loading is significantly different from this value, then a detector like BATSE 
will not be triggered. Dirty fireballs with Fq -C 300 will make long duration X- 
ray transients that will, in general, be too weak to trigger BATSE, and clean 
fireballs with Fq ^ 300 make brief high-energy 7-ray transients with insufficient 
fluence in the BATSE band to be detected. The implication is that there are many 
fireball transients that will be detected with more sensitive telescopes employing 
appropriate triggering properties and scanning strategies. 

No explanation has been given within the context of the colliding shell/internal 
shock model as to why GRB durations should range from a fraction of a seconds 
to hundreds of seconds. There seems to be no reason why intermittent or dealyed 
accretion of a massive ring of material around a collapsed star should not take place 
over long time scales, particularly given the unusual behavior that the accretion 
process must display if it is to produce the variability observed in GRB light curves. 
The observation of a single GRB that recurs after several hours, days, or months 
would falsify the external shock model. No convincing case of recurrence has been 
observed. 
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